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ON THE EXISTENCE OF PATHS BETWEEN POINTS IN 
HIGH LEVEL EXCURSION SETS OF GAUSSIAN 
RANDOM FIELDS 

By Robert J. Adler*'§ Elina Moldavskaya^^ and Gennady 

SAMORODNITSKY*'t 

The structure of Gaussian random fields over high levels is a 
well researched and well understood area, particularly if the field is 
smooth. However, the question as to whether or not two or more 
points which lie in an excursion set belong to the same connected 
component has constantly eluded analysis. We study this problem 
from the point of view of large deviations, finding the asymptotic 
probabilities that two such points are connected by a path laying 
within the excursion set, and so belong to the same component. In 
addition, we obtain a characterization and descriptions of the most 
likely paths, given that one exists. 

1. Introduction. Let X = (X(t), t G M. d ) be a real-valued sample 
continuous Gaussian random field. Given a level u, the excursion set of X 
over the level u is the random set 

(1.1) A u = {t£R d : X(t) > u) . 

Understanding the structure of the excursion sets of random fields is a math- 
ematical problem with many applications, and it has generated significant 
interest, with several recent books on the subject (e.g. [1] and [2]) and with 
considerable emphasis on the topology of these sets. One very natural ques- 
tion in this setting which has until now eluded solution but which we study 
in this paper is the following: given that two points in M. d belong to the excur- 
sion set, what is the probability that they belong to the same path connected 
component of the excursion set? Specifically, let a, b G R d , a / b. Recall 
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that a path in R d connecting a and b is a continuous map £ : [0, 1] — > R d 
with £(0) = a, £(1) = b. We denote the collection of all such paths by V(a, b) 
and are interested in the conditional probability 



P (3 £ G P(a, b) : X(£(v)) > u, for all < v < 1 X(a) > u, X(b) > uj . 

It is straightforward to check that we are considering measurable collections 
of outcomes, so this probability is well defined. 

Of course, the conditional probability above is a ratio of two probabili- 
ties, the denominator being no more than a bivariate Gaussian probability, 
which is well understood. Therefore, we will concentrate on the unconditional 
probability 

(1.2) *a,b(«) = P (3 i G V{a, b) : X(£(v)) > u, for all < v < 1) . 

If the random field is stationary, we may, without loss of generality, assume 
that b = 0, in which case we will use the simpler notation ^ a in (1.2). 

When the domain of a random field is restricted to a (compact) subset 
T cR d , the points a and b will be assumed to be in T, and the entire path 
in (1.2) will be required to lie in T as well. Nevertheless, we will use the same 
notation and also write 

^ a b (n) = P (3£ G V(a, b) : £(v ) G T and X(£(v)) > u, for all < v < 1) . 

Which of the two interpretations of is intended at any point will be 

clear from the context. 

We will study the logarithmic of the probability v I'a,b('w) for high levels u, 
i.e. as u — > 00. We start with a large deviations approach, which, as usual, 
will not only describe the probability but also give us insight into the highest 
probability configurations. This makes up Sections 3 and 4, which follow a 
brief technical Section 2 collecting some results on the reproducing kernel 
Hilbert space of a Gaussian process. Throughout we will treat the general 
and the stationary cases in parallel, but separately, since the stationary case 
is somewhat more transparent and more readily provides illustrative and 
illuminating special cases. In particular, we will look at a number of one- 
dimensional examples in Sections 5-7, where we can compute quite a lot. 
Even in this case the results are new and rather unexpected. We look at the 
multidimensional case in Section 8. While this section also contains some in- 
teresting and surprising examples, it turns out that typical examples involve 
non-convex optimization problems that we do not, at this stage, know how 
to solve in general. 
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2. Some technical preliminaries. In this section we introduce much 
of the notation we will use in the rest of the paper and recall certain impor- 
tant notions, concentrating in particular on the reproducing kernel Hilbert 
(RKHS) space of a Gaussian process. 

Our main reference for the RKHS is van der Vaart and van Zanten [9], 
and we use it selectively so as to prepare the background for using the large 
deviations theory of Deuschel and Stroock [3]. An alternative route would 
be to have followed the new notes by Lifshits [6]. 

We consider a real-valued centered continuous Gaussian random field X = 
(X(t), t G R d ). When needed (particularly, in the nonstationary case) we 
may restrict the domain of the random field to a compact subset T of M. d . 
We denote the covariance function of X by i?x(s,t) = cov(X(s), X(t)). 

As is customary, when the random field is stationary, we will use the single 
variable notation i?x(t) = -Rx(0,t) for the covariance function. In this case 
we denote the spectral measure of X by Fx, this being the symmetric, finite, 
Borel probability measure on M. d satisfying 



(2.1) i?x(t)= f e^Fx^x) 



t G 



If X is stationary, then this and local boundedness imply that 

Km ^ = 

||t||— >OC \\t\\ 

with probability 1, so that almost all the sample paths of X belong to the 
space 

Co = \ui = (oj(t), t G M d ), continuous, such that lim cj(t)/||t|| = >. 

I- ||t||— >oo J 

Equipped with the norm 

(2-2) Mo = sup — ^ , 

teR d 1 \\ z \\ 

Co becomes a separable Banach space, with dual space 

Cq = \ finite signed Borel measures n on R d with / ||t|| ||/z||(cit) < oo \. 

We view the stationary random field X as a Gaussian random element of 
Co, generating a Gaussian probability measure /^x on that space. 

In the absene of stationarity, we will usually consider a continuous Gaus- 
sian random field X = (X(t), t G T), for a compact set T C R d . In that 
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case we view the random field X as a Gaussian random element in the space 
C(T) of continuous functions on T, equipped with the supremum norm, thus 
generating a Gaussian probability measure fix ° n C(T). 

The reproducing kernel Hilbert space (henceforth RKHS) % of the Gaus- 
sian measure fix (or of the random field X) is a subspace of Co or C(T), 
depending on the parameter space of X, obtained as follows. In the gen- 
eral case we identify % with the closure C in the mean square norm of the 
space of finite linear combinations X^j=i a i^(^j) °f the values of the pro- 
cess, a,j G R, tj G M d (or T) for j = 1, . . . , k, k = 1, 2, . . . via the injection 
£ -»• C(T) given by 

(2.3) H ^w H = (E(X(t)H) , t € T) . 

When X is stationary, the RKHS % can also be identified with the sub- 
space, of functions with even real parts and odd imaginary parts, of the I? 
space of the spectral measure Fx in (2.1), via the injection L 2 {Fx) — > Co 
given by 

(2.4) /i -»• S(h) = (J e i(t ^ h(x) F(dx), telj . 

We denote by (-, •)% and || • ||% the inner product and the norm in RKHS %. 
Since both injections described above are isometric, we have the important 
equalities: 

(2.5) E{H 2 ) = \\wh\\h- 

In the stationary case, these can be written somewhat more informatively as 

(2-6) |H|| 2(Fx) = / ||M*)H 2 ^x(dx) = \\S{h)\\li- 

We shall use these equalities heavily in what follows. 

Note that for every s G M. d , the fixed s covariance function R s = R(-,s) is 
in %, and for every wjj G W, and t G lUij(t) = (wh, Rt)y.i meaning that 
the coordinate projections are continuous operations on the RKHS. This is 
also the reproducing property of the RKHS. Note also that the quadruple 
(C(T),H,w, nx) in general, or (Co,H, S, fix) hi the stationary case, is a 
Wiener quadruple in the sense of Section 3.4 in [3]. 

In the sequel we will use the notation M + {E) (resp. M^(E)) for the 
collection of all Borel finite (resp. probability) measures on a topological 
space E. 
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3. The basic large deviations result. We start with a large deviation 
result for the probability \I> a that exists a path between a and b wholly within 
a connected component of an excursion set. 

Theorem 3.1. . 

(i) Let X = (X(t), t G T) be a continuous Gaussian random field on a 
compact set T C Mr 1 . Then 

(3.1) lim -^log* aib («) = --C x (a,b), 

where 

C x (a, b) = inf {EH 2 : H G C, and, for some £ G 7>(a, b), 

(3.2) f (u) G T and > 1, < u < 1, } . 



(ii) Let X = (X(t), t G Mr) be a continuous stationary Gaussian random 
field, with covariance function satisfying 

(3.3) limsupiix(t) < 0. 

l|t|H°c 

Then 

(3.4) lim \]og* H {u) = -JCx(a), 
w/iere 

C x (a) = infj/ ||/i(x)|| 2 Fx(eZx) : /or some £ G V(0, a) 

,i(«i»),x) ^ Fx ( dx ) > i ; < u < 1 



/ 



PROOF. We start with putting our problem into the large deviation setup 
for Gaussian measures of [3]. We will use the language of Part (i) of the 
theorem, but the setup for Part (ii) is completely parallel. Observe that for 
u > 

*a,b(«) = P(« _1 X G A) , 
where A is the open subset of C(T) given by 

A = A^ h = G C{T) : 3 £ G V(a, b) such that w(£(u)) > 1, < v < l}. 
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Therefore, by Theorem 3.4.5 in [3] we conclude that 
(3.5) 

— inf I(oj) < lim inf — ^ log ^a b (u) < limsup — ^ log ^a b(u) < — infiYu;), 



for the rate function / which, by Theorem 3.4.12 of [3], can be written as 
(3.6) -- 



oo if uj £ %. 

for uj G C(T). Then (3.5) already proves the lower limit statement 

liminf \ logtf ab (u) > -^C x (a,b), 

valid for both parts of the theorem. Therefore, it remains to prove the match- 
ing upper limit. Here the argument is more involved in Part (ii) of the theo- 
rem, since non-compactness of the domain of the field requires us to rule out 
the possibility of increasingly long ranging paths. We present the argument 
in this case. The proof for Part (i) is similar, and easier (since we do not 
have to worry about paths which 'escape to infinity' as in the following). 

As is common with large deviation arguments, although we know that 
A = A° ^ A, this is not, per se, important. All that we need show is that 
the uj in the set difference A \ A do not contribute to the infimum on the far 
right of (3.5). 

We start by checking that 

(3.7) AC ( p| (1-5)aW H (l-6)A 1 ) 

\0<5<1 / \0<5<1 / 

(in the sense of the usual multiplication of a set of functions by a real num- 
ber), where A\ C Co is given by 

A 1 = S C : for every r > there is t G R d with ||t|| > r 

and a path f G V(0,t) such that uj(£(v)) > 1, < w < 1 j. 

To see this, let uj £ A, so that there is a sequence uj n G A, n = 1, 2, . . ., with 
uj n — > uj in Co- Suppose first that there is r > such that for a subsequence 
rik f oo, for each k = 1, 2, . . . there is a path € 7^(0, a) satisfying ||£fc(v)|| < 
r and WTi t (6W) > 1, < u < 1. Given < (5 < 1, choose A; so large that 

IK fe -w||c < 6/(1 + r). 



PATHS IN GAUSSIAN EXCURSIONS 



7 



Then for every t 6 K d with ||t|| < r we have |w nfc (t) — oj(t)\ < 5, so that 
<«>(&(«)) > 1 - * for < v < 1, and w € (1 - 5) A 

Alternatively, suppose that such r > does not exist. Then for every 
r > 0, for all but finitely many n, there is a path £ n 6 "P(0, a), going through 
a point t n with ||t n || = r, lying within the ball of radius r centered at the 
origin prior to hitting the point t n , and such that u} n (£ n (v)) > 1, < v < 1. 
Given r > and < (5 < 1, choose n outside of the above exceptional finite 
set, and so large that 

\\u n ~ u\\c <5/(l + r). 

As before, we conclude that there is a path connecting and t n such that the 
function oj takes values above 1 — S along this path. Therefore, co G (1 — S)A\, 
and so we have shown (3.7). 
Now note that since 



inf I(u) 

ue(l-5)A 



(1 



5f inf I(u) 

weA 



for any < S < 1, the rightmost inequality in (3.4), and so the result, will 
follow from (3.7) once we check that = oo for any w G ii, which we 

establish by showing that A\ D H = 0. 

Suppose that, to the contrary, there is a a; = S(h) G A for some h £ ~H. 
Fix an arbitrary e > 0. Assumption (3.3) guarantees the existence of a 
r £ > such that i?x(t) < e if ||t|| > r £ . By the definition of A, for every 
n = 1,2,... there is t„ with ||t n || = nr £ and a path £ connecting and t n 
such that u;(£(t;)) > 1, < v < 1. We can choose < v\ < . . . < v n < 1 
such that ||£(fj)|| = jr £ for j = 1, . . . , n. Then 



1 < 



if>(^)) = / \ l -jZ emV]U) I M x )^x(dx) 
n i=i Jli \ n H / 



< 



I^IIl 2 (f x ) 



i 2 (^x) 



However, 



Tl 

I e <«(«i).o 



7} 



i 2 (Fx) 



n— 1 n 

ni? x (0) + 2^ £ flx^*) 

Jl=l J2=jl+1 



< -fl x (0)+e, 
n 



8 



ADLER, MOLDAVSKAYA AND SAMORODNITSKY 



SO that 

IHli2(Fx)> ^x(0) + ^ 
Sending first n — > oo and then e — > we obtain ||/i||l2(^ x ) = oo, which is 
impossible. 

This contradiction proves the rightmost inequality in (3.4) and so we are 
done. □ 

Theorem 3.1 describes the logarithmic asymptotic of the path existence 
probability ^ a ,b m terms of a solution to an optimization problem in Hilbcrt 
space. The next result contains the dual version of this optimization problem, 
and relates ^ a ,b to the problem of finding a path of minimal capacity between 
a and b. 

Theorem 3.2. 

(i) Let X = (X(t), t G T) be a continuous Gaussian random field on a 
compact set T cR d . Then 

(3.8) lim ^log* a , b (u) 

= -^ x (a,b) 



sup min / / Rx(£(u) , £(«)) fi{du) fjt(dv) 

£eP(a,b)M6M+([0,l]) Jo JO 



(ii) Let X = (X(t), t G M d ) be a continuous stationary Gaussian random 
field, with covariance function satisfying (3.3). Then 

(3.9) lim ^log^ a (u) 

U— >OD U 



sup min / / i?x(C( n ) ~~ £.( v )) n{du) fi(dv) 

5eP(0,a)At6M+([0,l]) JO JO 



Note that the space M^QO, 1]) is weakly compact, and the covariance 
function i?x is continuous. Therefore, for a fixed path £, the function 

Jo Jo 

is weakly continuous on compacts. Hence, it achieves its infimum, and it is 
legitimate to write 'min' in (3.8) and in (3.9). 
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PROOF of Theorem 3.2. The proofs of the two parts are only notation- 
ally different, so we will suffice with a proof for Part (i) only. We use the 
Lagrange duality approach of Section 8.6 in [7]. Writing 

C x (a,b)= inf C x (a,b;0, 

?GP(a,b) 

where, for £ G V(&, b), 

(3.10)Cx(a,b;f) = inf {EH 2 : H G C and > 1, < v < 1} , 

we see that it is enough to prove that for every £ G P (a, b) , 



(3.11) C x (a,b;£) 



mm 



nin / / Rx(£(u),£(v)) fi(du) n(dv) 
+([o,i]) Jo Jo 



To this end, let Z = C([0, 1]. Then P = {z G Z : z(u) > 0, < u < 1} is a 
closed convex cone in Z. Its dual cone P® C Z* can be naturally identified 
with M+([0, 1]). Fix £ G P(a,b), and define G : £->Zhy 

G{H) = Gz(H) = (1 - w H {S(v)), 0<v<1). 
Then 67 is, clearly, a convex mapping. We can also write 

(3.12) (C x (a,b;0) 1/2 =mf{(£# 2 ) 1/2 : H € C, G(H) G — p| , 

and so our task is now to show that (3.12) implies (3.11). 

Suppose first that the feasible set in the optimization problem (3.10) is 
not empty. Then there is H G C such that G(H) belongs to the interior of 
the cone — P, so by Theorem 1, p224 of [7] we conclude that 
(3.13) 

(C x (a,b;£)) 1/2 = max inf (EH 2 ) 1/2 + f G(H)(v) n(dv) , 

and we may use 'max' instead of 'sup' because an optimal fi G M+([0, 1]) 
exists. For a fixed \i G M + ([0, 1]) with total mass ||//|| we let fi = fi/\\fJ.\\ G 
M+([0, 1]). Then 



(3.14) 



inf 

Hec 



(EH 2 ) 1/2 + f 1 G(H)(v) fi(dv) 
Jo 



= llull + inf 
Hec 



{EH 2 ) 1/2 -M [ w H {Z{v)) ftdv) 
Jo 
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lull + inf a 

a>0 



1 



Mil su P / w H {£(v)) p,(dv) 

H£C: EH 2 = l J 



-oo if > [sup HeC . EH 2 =l jQW H ((,(v))fi(dv)] 1 
\n\\ if < [sup HeC . EH 2 =l jQW H ((,(v))fi(dv)} _1 



Therefore, 



(C x (a,b;e)) 1/2 = 



inf sup / wh(^(v)) fi(dv) 

/J.&M+ ([0,1]) H£C:EH 2 =1 JO 



and (3.11) follows, since by the reproducing property of the RKHS, for every 
/iGM+([0,l]), 

sup f w H {£,{v)) n(dv) = sup [ (w H ,Rx(ti(v),-)) H »( dv ) 

H&C:EH 2 =1J0 HeC:EH 2 =lJ0 



sup 

w&H: ||io||-h = 1 



u\ / R x (€(v),-)[i(dv) 

/ H 

1 ,1 \ 1/2 

/ i2x(£(uU(v)M<*u)M<ftO 
o Jo 



In the last step we have used the fact that 



w„ = I' Rx(av),-)v(dv)eH, 
Jo 



so the supremum of the inner product is achieved at w = w^/Ww^W^i, and 



w,, 



[ [ R^(uU(v))fi(du)fJi(di 
Jo Jo 



1/2 



This establishes (3.11) for the case that the feasible set in (3.10) is not 
empty. We now turn to the case in which this set is, indeed, empty. This 
will complete the proof of the theorem. In this case (3.11) reduces to the 
statement 



(3.15) 



= min 

/i€M+([0,l])7o Jo 



Rx{Z(u),£(y)) fi (du) fi (dv) = 0. 



Suppose that, to the contrary, > 0. Let fig G M 1 + ([0, 1]) achieve the min- 
imum value in the integral defining I*. Consider the continuous real valued 
function 

W(u)= f fl x (£(«U(w)W<fr>)> 0<«<1. 
Jo 
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If this function never vanishes then, by continuity and compactness, it is 
bounded away from zero, so a sufficiently large in absolute value multiple of 
the random variable in C given by 

h= f 1 x(z(v))n>(dv) 

Jo 

is feasible for the optimization problem (3.10), contradicting the assumption 
that the set of feasible solutions is empty. 

Hence, there is no £ [0, 1] such that W(uo) = 0. For < e < 1 define a 
probability measure in ([0, 1]) by 

He = (1 - e)fi + eS U0 , 

where S a denotes the point mass at a. Note that 

1(e) = f 1 f 1 R^(uU(v)) fieidu) fi e (dv) 
Jo Jo 

= (l-e) 2 I* + 2e(l-e)W{u )+e 2 Rx(0) 
= (l-£) 2 /*+e 2 #x(0). 

Since was assumed to be positive, we see that 

d 



de 



Is) 



<0, 



e=0 



which contradicts the minimality of I*. This proves (3.15) and so the theo- 
rem. □ 

Observe that an alternative way of stating the result of Theorem 3.2 is 

-l 

(3.16) C x (a,b)= inf min 

£e7>(a,b) fi£M- 



in / / i?x(t, s) fi(dt) n(d& 



where M^(£) is the set of all probability measures in M. d supported by the 
path £ (strictly speaking, by the compact image of the interval [0, 1] under 
£). For a fixed path £ G V(a,b), the quantity 



(3.17) 



Cx(a,b;0 = 



mm 



m -Rx(t, s) n(dt) fi(ds) 

Jt Js 



is known as the capacity of the path £ with respect to the kernel i?x; see [4]- 
Therefore, on can treat the problem of solving (3.16) as one of finding a path 
between the points a and b of minimal capacity. 
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4. Fixed paths and measures of minimal energy. The dual for- 
mulation (3.4) of the optimization problem required to find the asymptotic 
of the path existence probability ^ aj b( u ) involves solving fixed path £ opti- 
mization problems (3.10) or (3.11). For a fixed path we have the following 
version of Theorems 3.1 and 3.2. 

Theorem 4.1. . 

(i) For a £ G V(a,b) let 

* a ,b(n;e) = P(x(£(v)) >u,0<v<iy 

Then 

(4.1) lim -llog* a)b («;0 = -^Cx(a,b;0. 

(ii) The primal problem (3.10) can be rewritten in the form 

(4.2) C x (a,b;£) = M^EH 2 : H G C, E[X(£(v))H] > 1, < v < l}. 

Further, if the feasible set in (4.2) is non-empty, then the infimum in (4.2) 
is achieved at a unique G C. 

(Hi) The set o//j£ Af-^QO, 1]) over which the minimum in the dual 
problem (3.11) is achieved is a weakly compact convex subset of ([0,1]). 
Furthermore, if the feasible set in (4.2) is non-empty, then, for every fx G W^, 

(4.3) ^({0 < v < 1 : E[X(Z(v))Hz] > l}) = 0. 

(iv) Suppose that the feasible set in (4.2) is non-empty. Then for every 
£>0, 

(4.4) P ( sup 

\0<u<l 

as u — > oo. Here 

(4.5) x t (v) = E[X(i(v))H^], 0<v< 1. 

The probability measures /j, G are called capacitary measures, or the 
measures of the minimal energy; see [4]. 

PROOF of Theorem 4.1. Part (i) of the theorem can be proved in the 
same way as Theorem 3.1. The fact that the primal formulations (3.10) and 
(4.2) are equivalent is an immediate consequence of the definition of wh- 



-X(Z(V))- X! :(V) 



> e 



X(£(v)) >u, 0<v<l] 
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Suppose now that the feasible set in (4.2) is non-empty, and let H n G C, 
n = 1, 2 ... be a sequence of feasible solutions such that EH% — > Cx(a, b; £). 
The weak compactness of the unit ball in C shows that this sequence has 
a subsequential weak limit with EH^ = Cx(a, b;£). Since the set of 
feasible solutions is weakly closed, is feasible. The uniqueness of the 
optimal solution to (4.2) follows from convexity of the norm. 

Convexity and weak compactness of the set follow from the nonnega- 
tive definiteness and continuity of i?x; see e -g- Remark 2, pl60, in [4]. The 
statement (4.3) is a part of the relation between the dual and primal optimal 
solutions; see Theorem 1, p224, in [7]. 

For Part (iv) of the theorem, note that by the Gaussian large deviation 
principle of Theorem 3.4.5 in [3], 

lim sup — ^ log P (x(^(v)) > u, < v < 1, sup 

u->-oo U \ Q<v<l 

(4.6) < -- inf^EH 2 , H G C : E[X(£(v))H] > 1, < v < 1, 

sup \E[X(t(v))H] - E[X(Z(v))Ht] \ > e\ . 

0<v<l J 

Therefore, the statement (4.4) will follow from Parts (i) and (ii) of the theo- 
rem once we prove that the infimum in (4.6) is strictly larger than Cx(a, b; £). 
Suppose that, to the contrary, the two infima are equal. By the weak com- 
pactness of the unit ball in C and the fact that the feasible set in (4.6) is 
weakly closed, this would imply existence of H* feasible for (4.6) such that 
EHl = EH^. Since is not feasible for (4.6), we know that i?* ^ H^. Since 
H* is feasible for (4.2), we have obtained a contradiction to the uniqueness 
of proved above. This completes the proof of the theorem. □ 

REMARK 4.2. Theorem 4.1 has the following important interpretation. 
Assuming that the feasible set in (4.2) is non-empty, Part (iv) of the theorem 
implies that the nonrandom function in (4.5) is the most likely choice 
for normalized sample path u~ 1 X(t;(v)), < v < 1 along £, given that 
{A" (£(«)) > u, < v < 1}. Part (iii) of the theorem implies that the values 
of the random field along the path £ have to (nearly) touch the level u at the 
points of the support of any measure of minimal energy. In other words, the 
sample path needs to be 'supported', or 'held', at the level u at the points 
of the support in order to achieve the highest probability of exceeding the 
high level u along the entire path £. We will see explicit examples of how 
this works in the following section, when we more closely investigate the one 
dimensional case. 



u 



X(S(v))-xt(v] 



> £ 
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The duality relation of the optimization problems (4.2) and (3.11) imme- 
diately provides upper and lower bounds on Cx(a, b;£) of the form 



(4.7) 



f 1 ^ Rx(£(uU(v))lA(du)lA(dv) 

Jo Jo 



<Cx( a ,h;0<EH 2 , 



for any fi G M{ l "([0, 1]) and any H 6 C feasible for (4.2). In particular, if 



(4.8) 



1 rl 



JO 



R x (Z(u)Z(v))v(du)v(dv) 



= EH 2 



for some \i and H as above, then [i 6 Wg, H = H^, and the common value 
in (4.8) is equal to Cx(a, b;£). 

Finding a measure of minimal energy, fj, £ W^, is, in general, a difficult 
problem. The following theorem includes a characterization of these mea- 
sures. 

THEOREM 4.3. Assume that the feasible set in (4.2) is non-empty. 
(i) For every fi G we have 



fl- c =C x (a,b;0 ^ X(Z(v))»(dv) 
Jo 



with probability 1. 

(ii) A probability measure G M-^QO, 1]) is a measure of minimal energy 
(i.e. fj, £ W^) if and only if 



(4.9) mm C R^{u)^{v)) fi(du) 

0<v<l Jq 



1 /•! 



JO 



Rx(£(ui),€(u2)) fJ-(du!) n(du 2 ) > 0. 



Note that Part (ii) of the theorem also says that the integral in the left 
hand side of (4.9) is equal to the double integral in its right hand side for 
//-almost every < v < 1. 

PROOF of Theorem 4.3. For Part (i), let G 14^. The calculations 
following the maximization problem (3.13) show that /i£ = Cx(a, b; ^) 1 / 2 /i 
is an optimal measure for that problem. It follows from Theorem 1, p224 in 
[7] that solves the minimization problem in (3.14), when any measure in 
M + ([0, 1]) optimal for (3.13) is used. Using the measure fi^ we see that 

H^a f 1 X(£(v)) n(dv) = aC x (a,b;0 1/2 f ' X(£(v)) fi(dv) 
Jo Jo 
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for some a > 0. Testing all random variables of the type 



H{ = b 



Jo 



f 1 X(£(v))»(dv), b>0, 



in (4.2) leads to the conclusion that 



i 



b 



min / 

0<«<i Jo 



The fact that b = Cx(a, b; £) follows now from the optimality of [i and the 
general properties of measures of minimal energy for bounded symmetric 
kernels; see e.g. Theorem 2.4 in [4]. 



We now prove Part (ii). Suppose first that fx satisfies (4.9), and define 
H e C by 



Jo 

where K(/j,) is the double integral in the right hand side of (4.9). Note that 
for any < v < 1, 



so that n and H satisfy the relation (4.8). Hence n 6 (and H = H^). 

In the opposite direction, if /x € W^, then the equality in (4.9) is a general 
property of measures of minimal energy for bounded symmetric kernels, as in 
the proof of Part (i). The fact that the equal terms in (4.9) are positive follows 
from the fact that the feasible set in (4.2) is non-empty, so Cx(a, b;£) < 
oo. □ 

REMARK 4.4. If the random field X is stationary, then the results of this 
section can be restated in the language used in Section 3 in the stationary 
case. In particular, the primal problem (4.2) becomes 





by (4.9). Therefore, H is feasible for (4.2). However, 




Cx(a;0 



= inf 



{/ ||/»(x)|| a F x (dx): / e 
{Jm. d jR d 



^)' x )/t(x)F x (dx) > 1, 0< v < 1 
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while the optimal solution of the primal problem in Part (i) of Theorem 4.3 
becomes 

/if(x) = C x (a;f) f e i{ ^ x) n(dv) for F x -almost all x € M. d , 
Jo 

for any fi G Wg. This relation can also be restated in terms on the measures 
supported by the (image of) path £ instead of the unit interval, as in (3.17). 
If fi is an optimal measure in (3.17), than we have 

(4.10) h ( (x) = C x (a; ^e i(t ' x) n(dt) for F x -almost all xGR* 1 . 

Note that the function in the right hand side of (4.10) is, up to a constant, the 
characteristic function of the measure [i. If the support of the spectral mea- 
sure Fx happens to be the entire space M. d , then the characteristic functions 
of all optimal measures in (3.17) are equal and, hence, the uniqueness of a 
characteristic function shows that, in this case (and as long as the feasible set 
in (4.2) is non-empty) there is exactly one probability measure /j, G (£) 
of minimal energy. 

REMARK 4.5. An immediate conclusion of Part (ii) of Theorem 4.3 and 
the assumed continuity of the covariance function is that the function 

v-> / Rx{£(u),Z(v))n(du),0<v<l, 
Jo 

is constant on the support of any measure ji 6 W^. This seems to indicate 
that the support of any measure of minimal energy may not be 'large'. In 
the examples below, however, this intuition holds only in some cases. 

5. The one-dimensional case. In this and the following two sections 
we specialize to the one-dimensional case d = 1. Let a < b. As before, we are 
interested the probability 

*a,b( u ) = P ( X ( t ) > u ' a <t < b ) ■ 

There is, essentially, a single path between a and b, and the results of the 
previous two sections immediately specialize to yield the following special 
case. (Note that condition (3.3) is superfluous in the one-dimensional non- 
stationary case.) 
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Theorem 5.1. Let X be a continuous Gaussian process on an interval 
including [a,b]. Then the limit 



-^C x (a,6)= lim -^log* a , 6 (u) 

2 U— >OD U Z 



exists, and 



C x (a,6) 

(5.1) = ini^EH 2 : H € £, E[X(a + (b - a)v)H] > 1, < v < l} 

(5.2) = min / / Rx(a + (b — a)u, a + (b — a)v)) [i(du)[j,(di 

neM+([o,i\)Jo Jo 

If the process X is stationary, an alternative expression for Cx(a) = Cx(0, a), 
a > 0, is given fry 

(5.3) 

r />oo poo ~\ 

C x (o) =infj/ ||/i(x)|| 2 F x (cfe) : / e™ a:c /i(x) F x (dx) > 1, < u < 1 \. 

The set W a ^ of n € M^ 1 " ([0, 1]) over which the minimum in (5.2) is achieved 
is a weakly compact convex subset of M± ([0,1]). The measures in W a ,b are 
characterized by the relation 

(5.4) Q m ^ 1 J R^.(a + (b — a)u, a + (b — a)v) fi(du) 

= / / Rx_(a + (b — a)ui,a + (b — a)u2) n(du\) p,(du2) ■ 
Jo Jo 

Suppose, further, that the problem (5.1) has a feasible solution. In this case 
the double integral in (5.4) is positive for any (i G W a ,b> an d the problem (5.1) 
has a unique optimal solution, H a jj. For each \x £ W a ^, 

H a)b = Cx(a,b) [ X (a + (b - a)v) fi(dv) 
Jo 

with probability 1. In the stationary case, the problem (5.3) has a unique 
optimal solution, h a . For each fi G W a = Wo,a 

h a (x) = Cx(a) / e iavx fx(dv) for F^-almost all -co < x < oo. 

JO 
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The conditional law on C[a, b] of the scaled process u~ 1 X restricted to the 
interval [a,b], given that X(t) > u, a < t < b, converges as u — >■ oo to the 
Dirac measure at 



(5.5) x a>b {t) = C x (a,b) R x (t, a + (6 - a)v) n(dv), a < t < b , 



Finally, if the process X is stationary, and the support of Fx. is the entire 
real line, then the set W a consists of a single probability measure, fj, a . 

Remark 5.2. Suppose that the process X is stationary. For /j, £ M{ l "([0, 1]) 
define fi = [ioT^ 1 with T : [0, 1] — > [0, 1] being the reflection map Tx = 1—x, 
< x < 1. If (i € W a , then jl satisfies conditions (5.4) because fi does, hence 
ft G W a as well. By convexity of W a , so does the symmetric (around x = 1/2) 
probability measure l/2(p, + ft). Therefore, W a always contains a symmet- 
ric measure. In particular, if W a is a singleton, then the unique measure of 
minimal energy is symmetric. 

In the remainder of this section we concentrate on the stationary case. 
We will investigate how the probability measure fj, a , the function h a and the 

limiting shape x a = Xo, a change as functions of a. This will help us under- 
stand the order of magnitude of the probability ^ a (u) for varying lengths a 
of the interval and, according to Part (iv) of Theorem 4.1, it will tell us the 
most likely shape the process X takes when it exceeds a high level u along 
the entire interval [0, a]. 

Our first result describes the situation occuring for some, but not all, 
stationary Gaussian processes on short intervals. 

PROPOSITION 5.3. Let X be a stationary continuous Gaussian process. 
Suppose that for some a > the following condition holds. 

(5.6) i2 x (t) + i? x (a -t)> i?x(0) + i?x(«) > for all < t < a. 

Then a measure in W a is given by 




and 



H ({0 < v < 1 : x a , b (a + (6 - a)v) > l}) = . 



(5.7) 



(1) A 1 x l x 



Furthermore, 



(5.8) 



Cx(a) = 



2 



^x(0) + i?x(«)' 
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_|_ e iax 

(5.9) h a {x) = — — — — -— for F-^-almost all — oo < x < oo, 

and 

(5 - 10) Xa{t)= i?x(0)+i?x(a) ' °-'- a - 

PROOF. Once we show that fj,^ £ W a , the rest of the statements will 
follow from Theorem 5.1. In order to prove (5.7) we need to check conditions 
(5.4). These follow immediately from (5.6) and the fact that 

jf 1 jf * Rx{a{ Ul - u 2 )) ^\d Ul ) ^\du2) = ^x(O) + l -R*{a) , 

while for < v < 1, 
-l 



R x {a(u - v)) n {1) (du) = -Rx(av) + -R*{a - av) . 



□ 



Remark 5.4. Note that a sufficient (but not necessary) condition for 
(5.6) is concavity of the covariance function i?x on the interval [0, a]. Indeed, 
for a concave covariance function the derivative exists apart from a countable 
set of points, and is monotone. Therefore, 

Rx{t) - Rx(0) = [ R'x(s)ds> I R' K (a-t + s)ds = Rx(a)-Rx(a-t). 
Jo Jo 

In particular, if the process X has a finite second spectral moment, then 
the second derivative of the covariance function exists, is continuous, and 
negative at zero (unless the covariance function is constant). Therefore, the 
derivative stays negative on an interval around the origin, hence the covari- 
ance function is concave on [0, a], and (5.6) holds, for a > small enough. 

On the other hand, apart from degenerate cases, the situation described in 
Proposition 5.3 cannot continue to hold for arbitrarily large a. For example, 
if the covariance function vanishes at infinity, then (5.6) fails for a large 
enough and t = a/2, say. 

In addition, a simple calculation shows that it always true that 

(5.11) lim log P (X(0) > u, X(a) > u) = -(i? x (0) + i?x(a)) _1 

u— >oo u 

Combining this with (5.8) shows that, in the scenario of Proposition 5.3, 
the probability that X exceeds a high level over an entire interval, and the 
probability that it does so only at the endpoints of the interval are, at a 
logarithmic scale, the same. 
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The plots of Figure 5.1 show the limiting shape x a for the stationary 
Gaussian process with covariance function Rx(t) = exp(— 1 2 /2), for a range 
of a for which Proposition 5.3 applies. In this case the largest such a is 
approximately equal to 2.2079. See Example 6.1 for more details. 



a=1.5 a=2.0 




0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 

t t 



Fig 5.1: Limiting shapes x a for the stationary Gaussian process with covari- 
ance function Rx(t) = exp(— 1 2 /2) when Proposition 5.3 applies. 

The plots of Figure 5.1 indicate that, as a approaches a critical value (ap- 
proximately 2.2079 in this case) the limiting curve x a 'attempts' to cross the 
level 1 at the midpoint of [0, a]. Equivalently, the normalized process u X 
attempts to drop below level 1 at that point and so, speaking heuristically, 
it has to be 'supported' at the midpoint t = a/2. The interpretation of The- 
orem 4.1 in Remark 4.2 calls for adding a mass to the measure /iW for the 
critical value of a at the midpoint of the interval. The next result shows that, 
in certain cases, this is indeed the optimal thing to do. 

PROPOSITION 5.5. Let X be a stationary continuous Gaussian process. 
Suppose that, for some a > 0, 



(5.12) 



Rx(P) + Rx(a) >2# x (a/2), 
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and let 

/r i n\ e i?x(0)+fix(«)-2fix(a/2) _ , 

(5 3) £a ~ 3i? x (0) + flx(a) - 4i? x (a/2) ( °' 1J " 

Suppose that for all < t < a/2, 

(5.14) i? x (i) + ^x(« - t) - i?x(0) - i?x(o) 

> e [#x(f) + #x(a - t) - i? x (0) - flx(o) - 2(i? x (a/2 - 1) - flx(a/2)) 
T/ien a measure in W a is given by 

(5.15) ^ i ^p^o + + £o * 1/2 . 
F-uri/iermore, 

3i?x(0) + i?x(a)-4i? x (a/2) 



(5.16) C x (a) 



(5.17) h a {x) = C x (a, 9 

/or Fx_-almost all — oo < x < oo, anc? 

1 - £a 



i?x(0) 2 + ^x(0)i? x (a) - 2i? x (a/2) 2 ' 
~1 -e a 



(l + e ia:E ) + e a e iaa;/2 



(5.18) z a (t) = C x (a) 
< t < a. 



i?x(i) + i?x(a - t)^+e a R x (\t - a/2| 



PROOF. The proof is identical to that of Proposition 5.3 once we observe 
that, under (5.12), // 2 ) is a legitimate probability measure. □ 

The plots of Figure 5.2 show the limiting shape x a for the stationary 
Gaussian process with covariance function Rx(t) = exp(— 1 2 /2), for a range 
of a for which Proposition 5.5 applies. In this case the range of a is, approx- 
imately, between 2.2079 and 3.9283. See Example 6.1 for more details. 

6. Specific covariance functions. In the previous section we saw some 
general results for one-dimensional processes, with some illustrative figures 
for what happens in the case of a Gaussian covariance function. In this sec- 
tion we look more carefully at this case, and also look at what can be said 
for an exponential covariance. 
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Fig 5.2: Limiting shapes for Rx(t) = exp(— t 2 /2) when Proposition 5.5 ap- 
plies (the top row). The left plot in the bottom row is a blowup of the 
right plot in the top row. The right plot in the bottom row shows how the 
constraints are violated soon after the upper critical value of a. 



EXAMPLE 6.1. Consider the centered stationary Gaussian process with 
the Gaussian covariance function 

(6.1) R{t) = e~ t2/2 , tel. 

For this process the spectral measure has a Gaussian spectral density which 
is of a full support in R. In particular, for every a > there is a unique 
(symmetric) measure of minimal energy. Furthermore, the second spectral 
moment is finite, so that, according to Remark 5.4, for a > sufficiently 
small this process satisfies the conditions of Proposition 5.3. To find the 
range of a for which this happens, note that conditions (5.6) become, in this 
case, 

(6.2) e~* 2 / 2 + e -( a -') 2 / 2 > 1 + e -« 2 / 2 , < t < a . 
Since the function 

g(t) = e- t2 / 2 + e- (a -^^,0<t<a, 
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is concave if < a < 2, and has a unique local minimum, at t = a/2, when 
a > 2, it is only necessary to check (6.2) at the midpoint t = a/2. At that 
point the condition becomes 

^(a) = 2e- a2 / 8 -l-e- a2 / 2 >0. 

The function ip crosses at a\ ~ 2.2079, which is the limit of the validity 
of the situation of Proposition 5.3 in this case. The plots of Figure 5.1 show 
the limiting shape x a for this process in the situation of Proposition 5.3. 

Somewhat longer (and numerical) calculations show that the conditions 
of Proposition 5.5 hold for the process with the covariance function (6.1) 
for an interval of values of a after the conditions of Proposition 5.3 break 
down. The conditions of Proposition 5.5 continue to hold until the second 
derivative at the midpoint t = a/2 of the limiting function in (5.18) becomes 
negative (so that the function takes values smaller 1 than in a neigbourhood 
of the midpoint). To find when this happens, we solve the equation 

(^x(i) + Rx(a ~ t))+e a R^(\t - a/2\) = 
at t = a/2. The resulting equation 

(l- £a )(a 2 /4-l)e- a2 / 8 - £a = 

has the solution a2 ~ 3.9283, which is the limit of the validity of the situation 
of Proposition 5.5 in this case. The plots of Figure 5.2 shed some light on the 
above discussion. This discussion indicates, and calculations confirm, that, 
in the next regime, the mass in the middle for the optimal measure splits 
into two parts that start to move away from the center. Heuristically, this 
is needed "to support" the trajectory that, otherwise, would "dip" below 1 
outside of the midpoint. 

These calculations rapidly become complicated. They seem to indicate 
that the next regime continues to hold until around 03 « 5.4508. In this 
regime the optimal measure takes the form 

(6-3) M (3) = + l ~=^5i + §<5 V 2-<*a + f <W> 

where d a is the distance of two internal masses from the midpoint. When 
a = 4.5, £4.5 = 0.36632 and 04.5 = 0.12285, so that the internal atoms are 
at 0.37715 and 0.62285, and the rest of the support is concentrated at the 
endpoints of the interval with probabilities 0.31684. Figure 6.1 shows the 
limiting shape X4.5. 
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a=4.5 




12 3 4 



Fig 6.1: The limiting shape in the case a = 4.5 for i?x(i) = exp(— i 2 /2) 

It would be nice to understand all regimes, but we do not yet know how 
to find a general structure. On the other hand, Section 7 gives asymptotic 
results for a — > oo. 

Finally, Figure 6.2 shows the growth of the exponent Cx(a) with a for as 
long as either Proposition 5.3 or Proposition 5.5 applies. 



8 




12 3 



Fig 6.2: The exponent Cx(a) as a function of a for Rx(t) = exp(— 1 2 /2). 

The next example shows a situation very different from that of Example 
6.1. 

EXAMPLE 6.2. Consider an Ornstein-Uhlenbeck process, i.e. a centered 
stationary Gaussian process with the covariance function 

(6.4) R(t) = e -1 * 1 , tel. 
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For this process the spectral measure has a Cauchy spectral density, so it 
is also of full support in R. Therefore, for every a > there is a unique 
(symmetric) measure of minimal energy. In this case, however, even the first 
spectral moment is infinite. The covariance function is actually convex on the 
positive half-line so, in particular, the conditions of Proposition 5.3 fail for 
all a > 0. In fact, it is elementary to check that for the probability measure 

(6-5) n = -^-:5o + — + -J- A 

a + 2 a + 2 a + 2 

(where A is Lebesgue measure on (0, 1)), the integrals 



L 



i 

Rx_{a{u — v)) fi(du), < v < 1 



have a constant value, equal to 2/(a + 2). Therefore, the measure fi in (6.5) 
is the measure of minimal energy, and Cx(a) = (a + 2)/2 for all a > 0. 

By Theorem 5.1 we conclude that the limiting function x a is equal to 1 
almost everywhere in [0, a] with respect to Lebesgue measure. Since x a is 
continuous, it is identically equal to 1 on [0, a]. 

Examples 6.1 and 6.2 demonstrate a number of the ways a stationary 
Gaussian process 'prefers', in the large deviations sense, to stay above a high 
level over an interval. The process of Example 6.1 with covariance function 
(6.1) is smooth; the likeliest way it stays over a level is to force it to be 
'slightly' over that level at a properly chosen finite set of time points; after 
that it is 'held' over the level at the rest of the interval [0, a] by the corre- 
lations of the process. The optimal configuration of the finite set of points 
depends on length of the interval [0, a], and it appears to undergo phase 
transitions at certain critical interval lengths. The complete picture of this 
'dynamical system' of finite sets remains unclear. On the other hand, the 
Ornstein-Uhlenbeck process of Example 6.2 is continuous, but not smooth. 
In fact, it locally behaves like a Brownian motion. Therefore, 'holding' it 
'slightly' above a level at a discrete point does not help, since it 'wants' im- 
mediately to go below that level. This explains the nature of the optimal 
measure jjl in (6.5), and this nature stays the same no matter how short or 
long the interval [0, a] is. In particular, phase transitions do not happen for 
this process. 

It remains to be investigated whether other types of behavior are possi- 
ble, and under what exact conditions on the Gaussian process each type of 
behavior occurs. It is also likely that minimal energy measures in W a carry 
additional information, describing how 'slightly' over the level u a Gaussian 
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process is most likely to be given that it is over that level along the interval. 
The exact nature of this information also remains to be investigated. 

7. Asymptotics for long intervals. In this section we investigate the 
asymptotics of the exponent Cx(a) for large a. We start with a result showing 
that, for certain short memory stationary Gaussian processes, the exponent 
Cx(a) grows linearly with a over long intervals. Furthermore, the energy of 
the uniform distribution A on [0, 1] becomes, asymptotically, minimal. 

Theorem 7.1. Let X be a stationary continuous Gaussian process. As- 
sume that i?x is positive, and satisfies the following condition: 

roo 

(7.1) / R(t)dt <oo. 

Jo 

Then, with X denoting the uniform probability measure on [0, 1], 

(7.2) lim 1 C x (a) = ( lim a I [ R^(a(u - v)) X(du)X(di 

a— >oo a ya^-oo Jo Jo 



2JZ°R(t)df 

PROOF. By Theorem 5.1, the statement of the the present theorem is 
equivalent to the following pair of claims. 

rl rl f'OO 

(7.3) lima/ / R x (a(u - v)) X(du)X(dv) = 2 R(t) dt 
a ^°° Jo Jo Jo 

and 

rl rl poo 

(7.4) lim info min / / R^(a(u- v)) n{du)n{dv) > 2 / R(t)dt. 

a ^°° M eM+([o,i]) Jo Jo Jo 

Since 

»i r i 

Rx_(a(u - v)) X(du)X(dv) 



o Jo 



2 rl rav ra(l—v) 

- / R(t)dt+ / R(t)dt 
a Jo Jo Jo 



dv , 



(7.3) immediately follows from (7.1) and bounded convergence theorem. 
Therefore, it only remains to prove (7.4). Suppose that, to the contrary, 

(7.4) fails, and choose a sequence a n — > oo such that 

rl rl roo 

lim a n min / / R^_(a n (u — v)) fi(du) ix(dv) < 2 / R{t)dt. 
n ^°° M6M+([0,1]) Jo Jo Jo 
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For each n choose a symmetric fi an 6 W an , so that 



(7.5) lim a„ / / R^(a n (u - v)) fi an (du)n an (dv) < 2 R(t) dt . 
n ^°° Jo Jo Jo 

We claim that for every 7 > 0, 

(7.6) lim Man ([0,7a; 1 ]) =0. 



Indeed, by the positivity of i?x, for any 7 > 
•1 r l 



f f Rx{a n (u-v)) n an (du)[j, an (dv) 
Jo JO 

pa-n 1 pan 1 

> / / Rx(a n (u- v)) fi an (du)fi an (dv) 
Jo Jo 



2 

.-11 ' 



so that (7.5) necessitates (7.6). Next, define a sequence of signed measures 
on [0, 1] by p, n = fi an — X. Note that 

(7.7) An([0,l]) = for each n. 

By the nonnegative definiteness of i?x, 

[ [ i?x (a n {u - v)) Ha n (du)n an (dv) 
Jo Jo 

= / / Rx{a(u - v)) X(du)X(dv) + / / R^(a n (u - v)) fi n (du)fi n (dv) 
Jo Jo Jo Jo 

+2/ / Rx(a n (u - v)) X(du)fi n (dv) 
Jo Jo 

I / Rx{a(u - v)) X(du)X(dv) + 2 / / R^(a n (u - v)) X(du)p, n (dv). 
Jo Jo Jo 



We will show that 

-1 r-i 



(7.8) lim a n / / i2 x (a„(u - v)) X(du)p, n (dv) = . 

Together with (7.3) this will provide the necessary contradiction to (7.5). 
Let 7 > 0. Write the integral in (7.8) as 



/•i-7Cn 1 r rl 

/ / Rx(a n (u ~ v)) du 

J-ya^ 1 Uo 



P>n(dv) 
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rja-n 1 r /•! 

+2 / / Rx(a n (u - v)) du p, n (dv) 

Jo Uo 

± W + 2 J« • 



Observe that 



j; 



(2)1 



\ flo, n ra„v ra n {i—v) 

— / Rx{t)dt+ / 22 x (t) 

0"a JO Jo Jo 

2j °°flx(t)rft|,. ,,/ rn _n\ 

so that by (7.6) we obtain 

(7.9) lim OnJj® = 



dt 



fln{dv) 



for every 7 > 0. Next, we write 



1 /•!— 7a„ fa„» /•a„(l— u) 

4 1} = — / / Rx(t)dt+ / fl x (t) 

«n J^a" 1 [JO JO 

2j °°i?x(t)dt , / x 

= — y ^n ( [7«„ . 1 - 7«n J ) 

fl ra V / 

X /"l-T^n 1 [ /"oo roc 

/ / Rx{t)dt+ / R x (t) 

On J'ya n 1 J a n v J a n (l— v) 

A /( ii) _ ,(12) 

It follows from (7.7) that 

(7.10) \a n 4 u ^\=2^ /? x (i)^|An([0,7On 1 ])+An([l-7an 1 >l 
as n — >■ 00, by (7.6). Finally, 



fit 



fln(dv) 



dt 



P>n(dv) 



(7.11) 



/>oo 

|«nJ n 12) |<4/ i?x(t)A, 



and we obtain by (7.9), (7.10) and (7.11) that 



lim sup a n 

n^oo 



r 1 />1 /•oo 

/ / R*(a n (u-v))\(du)ji n (dv) <4/ fl x (t)dt. 

JO JO J7 



Letting 7 — >■ 00 proves (7.8) and, hence, completes the proof of the theorem. 

□ 
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The next theorem is the counterpart of Theorem 7.1 for certain long mem- 
ory stationary Gaussian processes. In this case, the uniform distribution on 
[0, 1] is no longer, asymptotically, optimal. We will assume that the covari- 
ance function of the process is regularly varying at infinity: 

(7.12) R*(t) = ^ , < /3 < 1, 

where L is slowly varying at infinity. Before stating the theorem, we introduce 
new notation. 

Consider the minimization problem 

(7.13) min C C ^(dv) < /3 < 1. 
ixeM+ ([o,i]) Jo Jo \u-v\p 

This is a minimization problem of the same nature as in (5.2) with a = 0, b = 
1, and the covariance function i?x replaced by the Riesz kernel Kp(u,v) = 
\u — v\~P, u, v £ [0, 1]. The general theory of energy of measures in [4] applies 
to the Riesz kernel. In particular, the minimum in (7.13) is well defined, is 
finite and positive. Let Wp be the set of measures in M+ ([0, 1]) of minimal 
energy with respect to the Riesz kernel. Note that the uniform measure 
A Wp since it does not satisfy the optimality conditions in Theorem 2.4 in 
[4]- 

THEOREM 7.2. Let X be a continuous stationary Gaussian process. As- 
sume that Rx is positive, and satisfies assumption (7.12) of regular variation. 
Then for any fip G Wp, 

(7.14) Mm fixMCx(a) = (ff ) " . 

a^oo VJ Jo \U-V\P J 

PROOF. Suppose first that there is a sequence a n f oo such that 

(7.15) lim — — - — - min / / Rji_(a n (u — v))fi(du)fi(dv) 

n-Kx. i? X (On) t*eM+ ([0,1]) Jo Jo 

Z 1 f 1 Hp{du)np(dv) 



o \u-v\P 



For each n choose \x n G W an , let n& f 00 be a subsequence such that fj, nk ft, 
weakly as k — > oo for some ft G M 1 + ([0, 1]). By Fatou's lemma and the regular 
variation of i?x> 

n , millf P In \ l l R ^i a n k {u-V)) fl nk (du)n nk (dv) 
fc^oo Kyi{a nk ) Jo Jo 
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>-ff 

Jo Jo 



p,(du) jl(dv) 



\u — v\P 
1 f 1 /j,p(du)fip(dv) 
lo Jo \u-v\f> 



since \ip has the smallest energy with respect to the Riesz kernel. This con- 
tradicts (7.15), thus proving that 

(7.16) liminf — — — min / / Rx_(a(u — v)) fi(du)fi(dv) 

a^oo i? x (a) peM+([0,l])Jo Jo 

1 



> 



// 

Jo Jo 



\u — v\P 

In order to finish the proof, we need to establish a matching upper limit 
bound. 

To this end, let 9 > be a small number. We define a probability measure 
ve G M^([0, 1]) by convolving up with the uniform distribution on [0, 9] and 
rescaling the resulting convolution back to the unit interval. More explicitly, 
if X and U are independent random variables, whose laws are jip and A, 
respectively, then ug is the law of (X + 9U)/(1 + 9). Note that 

(7.17) v e < A and ^ < a.e. on [0, 1]. 

aA 9 

Given < e < 1 — f3, by Potter's bounds (see e.g. Proposition 0.8 in [8]), 
there is to > sufficiently large to ensure 

for all t > to and x > 1. We have 

min / / Rx(a{u — v)) fi(du)fi(dv) 
fieM+ ([o,i]) Jo Jo 

< / Rx(a(u - v)) v e (du)v g (dv) 
Jo Jo 

f / l(|it — v I < to/a^Rx(a{u — v)) ug(du)i/g(dv) 
o Jo 

/ / l(\u — v\ > to/a)Rx(a{u — v)) ug(du)ug(dv) 
Jo Jo 



= h(a) + I 2 (a) . 
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By the definition of to, 

l(\u — v\ > to/a) 



Rx(a(u 



-Rx(a) ~ 1 + £ 
so that by the dominated convergence theorem we have 

1 vg(du)vg(dv) 
\u — v\@ 

o{Rx(a)) 



u-v\-^ 



On the other hand, by (7.17), 

2i 1 + 



as a 



/i(a) < Rx{0)'- 

oo. We conclude that 
1 



6 



1 /•! 



lim sup — — — min 
a^oo itx(a) M eM+([o,i])7o io 

'* 1 Ug(du)ug(dv 



Rx(a(u — v)) fi(du)/j,(dv) 



< 



II 

Jo Jo 



\u — v\P 



Once we show that 



(7.18) 



lim 

e-s>o 



1 rl 



•/ 



Vg{du)Vg{dv) 



\u — v\ 



1 rl 



JO 



Hp(du)np(dv) 



\u — v\ 



we will have established an upper bound matching (7.16). This will complete 
the proof of the theorem. Recall that (7.18) is equivalent to 



lim E 



X 1 - X 2 + 0(C/i - U 2 ) 



-a 




= E 


X\ — X 2 



where Xi, X 2 , U±, U 2 are independent random variables, X\ and X 2 with the 
law fip, while U\ and U 2 are uniformly distributed on [0, 1]. This, however, 
follows by the dominated convergence theorem and the following fact, that 
can be checked by elementary calculations: there is rp G (0, oo) such that for 
any < b < 1 and < 6 < 1, 



E 



b + 0(L/i - U 2 ) 



< rpb- 



□ 



REMARK 7.3. It follows from Proposition A. 3 in [5] that the energy of 
the measure \x$ with respect to the Riesz kernel cannot be smaller that one 
half of the energy of the uniform measure. Hence, 

lim i?x(a)Cx(a) G ((1 - 0)(2 - /3)/2, (1 - /3)(2 - /?)) . 
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8. The multi-dimensional case. Our understanding of the one-dimen- 
sional case described in the previous three sections, while incomplete, is nev- 
ertheless quite significant. In contrast, there is much less we can say about 
the multivariate problem of Section 3. The problem lies, in part, in the non- 
convexity of the feasible set in (3.1) which leads, in turn, to the 'max-min' 
problem in Theorem 3.2. 

The following proposition is a multivariate version of Proposition 5.6. Note 
that stationarity of the random field is not required. 

PROPOSITION 8.1. Let X = (X(t), t e T) be a continuous Gaussian 
random field on a compact set T C M. d , and suppose that a, b are in T. 
Suppose that there is path £o in T connecting a and b such that 
(8.1) 

p < ,mUb /fMKW fl x (a, a) + 2i? x (a, b) + fl x (b, b) 
i?x(a,£o(w)j +R x (£ (u),b) > > 

for all < u < 1. Then the supremum in (3.8) is achieved on the path £o 
and 

(8.2) C x (a,b) 



i?x(a, a) + 2i? x (a, b) + fl x (b, b) 



Remark 8.2. Using u = and u = 1 in (8.1) shows that conditions 
of Proposition 8.1 cannot be satisfied unless i?x(a, a) = i?x(b,b). Corre- 
spondingly, we can restate (8.2) as 

Cx(a,b) = 



i?x(a,a) + J R x (a,b) ' 

Recall (5.11), which shows that this implies the logarithmic equivalence of the 
probabilities of X being above the level u along a curve or at its endpoints. 

PROOF of Proposition 8.1. Consider the fixed path £ - The assump- 
tion (8.1) shows that the measure 

1 r 1 r 
IM) = 7,°a + 2° 6 

satisfies conditions (4.9) and, hence, is in W^ by Theorem 4.3. Therefore, 
min / / i?x(6W,6W) n{du)n{dv) 

(i£M+([0,l]) JO JO 

-1 r l 

#x(£o no(du) no(dv) 

o Jo 
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_ fix (a, a) + 2fi x (a, b) + i? x (b, b) 
4 

On the other hand, for any other path in T connecting a and b, 

min I [ Rx(£(u),£(y)) n(du)n(dv) 
neM+([o,i]) Jo Jo 

< [ I R*((i{u)^{v))no{du)iio(dv) 
Jo Jo 

_ fix(a,a) + 2fi x (a,b) + fi x (b,b) 

4 

Therefore, the supremum in (3.8) is achieved on the path £o> an d (8-2) follows 
by Theorem 3.2. □ 

Even for the most common Gaussian random fields, the assumptions of 
Proposition 8.1 may be satisfied on some path but not on the straight line 
connecting the two points. In that case, the straight line, clearly, fails to be 
optimal. 

EXAMPLE 8.3. Consider a Brownian sheet in d > 2 dimensions. This is 
a continuous centered Gaussian random field X on [0, oo)^ with covariance 
function 

d 

-Rx(s,t) = JJmin(sj,tj), s,t G [0, oo) d . 

3=1 

We restrict the random field to the hypercube T = [0, d] d , and let 

a= (l,2,...,d-l,d), b = (d,l,2, 

It is elementary to check that the path 

f (1 + d(d - l)u, 2, ...,d- 1, d) for < u < \ 

io{u) = \ (d,l,...,j-2,2j-l-du,j + l,...,d) forl^< u <i, 
{ j = 2,...,d 

satisfies (8.1) and, hence, the supremum in (3.8) is achieved on that path. 
Therefore, by Proposition 8.1, 

C x (a,b) = 



d! + (d-l)! 



On the other hand, if we consider the straight line connecting the points a 
and b, 

f (u) = (l + (d - l)u, 2 - u, 3 - u, . . . , d - u) , 0<u<l, 
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then the sum in the right hand side of (8.1) becomes 
d 

L{u) = Y[(j -u) + (l + (d- l)u)(d - 1)!, < u < 1 . 

i=2 

The function L achieves the value d\ + (d — 1)! at the endpoints u = and 
u = 1, and is strictly convex if d > 3. Therefore, it takes values strictly 
smaller than d\ + (d— 1)1 over < u < 1. That is, (8.1) fails, and the straight 
line is not optimal. If d = 2, however, then L is a constant function, condition 
(8.1) holds over the the straight line path, and the straight line is optimal. 

We also note that, if d = 1, then the Brownian sheet becomes the Brow- 
nian motion in one dimension. In that case it is, clearly, impossible to find 
two positive points a < b in which the process has the same variance, so 
Proposition 8.1 does not apply. In this case, however, we are in the situation 
of Theorem 5.1, so if < a < b < oo, then the measure (i = 5 a satisfies (5.4) 
and, hence, is optimal. 

The above example notwithstanding, under certain assumptions on the 
random field, the straight line path between two points turns out to be 
optimal for the optimization problem (3.8). The next result describes one 
such situation. 

Recall that a random field on M. d is isotropic if its law is invariant under 
rigid motions of the parameter space. A centered Gaussian random field X is 
isotropic if and only if its covariance function is a function of the Euclidian 
distance between two points. With the usual abuse of notation we will write 
i?x(a,b) = J R x (||b-a||), a,beR d . 

PROPOSITION 8.4. Let X be a continuous centered isotropic Gaussian 
random field, such that the covariance function i?x is nonincreasing. Then 
for any a, b € M. d , the straight path connecting the points a and b is optimal 
for the optimization problem (3.8). 

PROOF. We may and will assume, without loss of generality, that a = 
(a, 0, ... , 0) and b = for some a > 0. We start with showing that the 
supremum over £ G V(0, a) is achieved over paths in 

V { = |^ : [0,1] {(a:,0, ...,0) : x > 0}, continuous, £(0) = 0, £(1) = a} . 

To this end, it is enough to show that for each £ G "P(0,a) there is £ G V\ 
such that 

(8.3) min C #x(||£H - f (v)||) n(du) fi(dv) 

MeM+([o,i]) Jo Jo 
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r-l fl 



< min 



ain / / Rx(\\i(u) - £(v)\\) v(du) v(dv) 
+([o,i]) Jo Jo 



To see this, define for £ G P^a) 

|(iO = (||£(u)||,0,...,0), 0<n<l. 

Clearly, £ eVi, and (8.3) follows by the monotonicity of i?x and the triangle 
inequality 

|||f(u)||- < lieW -t(v)\\. 

Next, any £ G T 3 ; is of the form 

(8.4) e(«) = M«),o,...,o), 0<«<1, 

with </9 : [0, 1] — )■ [0, oo) a continuous function, satisfying <^(0) = 0, (p(l) = a. 
Defining <p(u) = min(<£>(u), a), < u < 1, and 

|(u) = (0(u),O,...,O), 0<u<l, 

we see that the supremum over paths in is, actually, achieved over paths 
whose image is exactly the interval [0, a]. Finally, for any path £ G of the 
latter type, given in the form (8.4), define 

r(v) = inf{u G [0, 1] : <p(u) = av}, < v < 1 . 

Then r is a measurable map from [0, 1] to itself, so for any fi G ([0, 1]), 
we can define fii G M+([0, 1]) by /ii = fi o r _1 . Then 

l /-l 

-Rx(|</>(^) - 9?( w )l) Mi(^) V>i{dv) 
l ,i 

/ fl x (|p(r(u)) - <p(r(v))\) »(du) ^(dv) 

Jo 

1 ,1 

/ ^?x(a|"U — v|) ^{du) n{dv). 
io Jo 

Therefore, 

1 % x (||^-^)llM^)M^) 

< min / / Ry^(a\u — v\) fj,{du) fi(dv), 
M eM+([o,i])7o 

and the statement of the proposition follows. □ 



o JO 



mm 

AteM+([o,i]) Jo Jo 
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